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1.  Introduction 

1.A  Background: 

Each  year  about  186,000  women  contract  breast  cancer,  of  whom 
46,000  will  die  of  it.  Breast  cancer  is  the  leading  cause  of  cancer 
deaths  amongst  American  women  aged  35  to  50  and  also  has  the 
highest  cancer  mortality  rate  for  younger  women.  Mammography 
is  most  widely  used  modality  for  early  detection  of  breast  cancer. 

The  detection  of  tiny  carcinomas  and  low  contrast  tissue  is  essen¬ 
tial  for  mammography  to  reduce  the  mortality  rate1,  2 .  It  has  been 
shown  that  about  20%  of  detectable  cancers  are  overlooked  or  mis¬ 
diagnosed  on  first  detection.  We  obviously  need  a  major  improve¬ 
ment  in  the  quality  of  mammograghic  imaging.  Using  an  innova¬ 
tive  new  technology  x-ray  optics  to  improve  system  resolution,  en¬ 
hance  contrast  and  reduce  dosage  will  increase  the  effectiveness  of  this  proven  screening  modality,  with 
a  direct  impact  on  mortality. 

1.B  Benefits  of  digital  mammography 

Conventional  film  screen  mammography  suffers  from  limited  dynamic 
range  and  film  granularity,  which  can  reduce  the  sensitivity  to  detect  micro¬ 
calcifications.  The  primary  theoretical  limitations  of  mammography  are  the 
system  resolution,  which  determines  the  minimum  size  of  the  detectable 
malignancy,  and  the  need  to  expose  the  patient  to  ionizing  radiation.  Using 
innovative  new  technology  —  polycapillary  optics  —  to  improve  system 
resolution  and  reduce  required  dose  will  increase  the  effectiveness  of  this 
proven  screening  modality,  with  a  direct  and  immediate  impact  on  mortality. 
By  matching  polycapillary  optics  to  digital  detectors,  we  can  develop  direct 
digital  mammographic  systems  which  can  avail  themselves  of  the  advan¬ 
tages  of  digital  processing,  including  improved  image  contrast  and  resolu¬ 
tion  at  reduced  radiation  dose.  In  practice,  mammographic  imaging  is  often 
limited  by  quality  assurance  issues,  which  can  also  be  favorably  addressed 
by  digital  processing.  Digital  detection  can  provide  high  dynamic  range 
which,  in  addition  to  improving  contrast,  greatly  increases  the  tolerance  of  the  final  image  to  under  or 
over  exposure.  Digital  images  can  be  enhanced  and  are  amenable  to  computer-aided  diagnosis.  Finally, 
digital  images  can  also  be  transported  quickly  for  skilled  consultations. 

1.  C  What  capillary  x-ray  optics  can  do 

Capillary  optics,  consisting  of  arrays  of  hollow  glass  tubes,  is  a  relatively  new  technology  for  controlling 
X-ray  beams.  X-rays  incident  on  the  interior  of  the  glass  tubes  at  small  angles  are  guided  down  the  tubes 
by  total  external  reflection.  Arrays  of  curved/tapered  capillaries  can  be  used  to  focus,  collimate  and  filter 
X-ray  radiation3,4,5,6,7.  Such  arrays  are  manufactured  by  stringing  hollow  glass  polycapillaries  through 
metal  grids  (as  in  Figure  1)  or  as  a  monolithic  optic. 


Figure  2  Cross  section 
of  polycapillary  fiber. 
0.5  mm  in  diameter. 
SEM  photo. 
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Figure  l  Complete  optic  con¬ 
structed  by  configuring  fibers 
through  metal  grids. 
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The  use  of  Kuinakhov  capillary  optics  in  place  of  conventional  scatter  reduction  grids  in  a 
mammographic  system  has  significant  potential  to  provide  improved  resolution,  increased  contrast  en¬ 
hancement  and  reduced  dose  in  mammographic  imaging.  The  optics  can  also  be  used  to  mate  the 
radiographic  image  with  a  digital  detector  by  appropriate  choice  of  demagnification  and  separation  to 
discrete  chips.  In  addition,  a  pre-patient  optic  could  be  employed  to  increase  the  available  intensity  in  a 
fan  beam  relative  to  simple  slot  collimating. 

The  fibers  in  this  research  varied  from  0.3  mm  to  4  mm  in  outer  diameter,  with  channel  sizes  (the 
diameters  of  the  holes  within  the  fibers)  ranging  from  4  to  22  pm,  as  shown  in  Table  1  in  section  2.E. 
Open  area  in  the  table  refers  to  the  fraction  of  the  front  face  of  the  fiber  which  is  open  channels,  as 
opposed  to  glass  walls.  No  fiber  can  have  a  transmission  that  is  greater  than  the  open  area. 

X-ray  beams  can  be  bent,  focused  or  collimated  by  a  carefully  curved  polycapillary  array.  X-ray  photon 
energies  are  much  larger  than  the  plasma  frequencies  of  glasses,  which  are  tens  of  electron  volts.  In  this 
regime,  the  real  part  of  the  index  of  refraction  of  glass  can  be  simply  approximated  by 


where  n  is  the  index  of  refraction,  e  is  the  dielectric  constant  of  the  glass,  to  is  the  photon  frequency,  top 
is  the  plasma  frequency  of  the  material,  and  so  is  the  dielectric  constant  for  vacuum. 

The  plasma  frequency  cop  for  glass  is  small  compared  to  the  photon  frequency,  so  that  n  is  slightly  less 
than  unity.  Therefore,  x  rays  traveling  in  vacuum  or  air  can  be  totally  externally  reflected  from  smooth 
glass  surfaces.  The  maximum  angle  between  the  x-ray  beam  and  glass  surface  for  total  reflection,  0C,  is 
approximately 

0  c  =  mrad,  (2) 

where  E  is  the  photon  energy  in  keV.  The  higher  the  energy  of  photon,  the  smaller  is  the  critical  angle. 
For  this  reason,  profile  errors  that  change  the  angle  of  incidence  are  more  detrimental  at  high  energies. 


1.D  Numerical  simulations 

To  evaluate  the  experimental  performance  of  capillaries  and  optics,  and  to  design  capillary  optics,  it  is 
necessary  to  be  able  to  predict  theoretical  behavior  for  complex  geometries.  Extensive  modeling  pro¬ 
grams  describing  the  propagation  of  x-rays  along  capillaries  and  optics  are  being  developed.  These 
modeling  programs  provide  essential  information  on  the  transmission  efficiency  and  divergence  of  cap¬ 
illary  optics,  and  are  being  used  to  simulate  the  performance  of  the  optics  by  obtaining  the  following 
outputs:  transmission  efficiency,  exit  divergence,  output  uniformity,  and  transmission  as  a  function  of 
entrance  angle. 


The  models  for  these  simulation  programs  are  based  on  Monte  Carlo  simulations  of  simple  geometrical 
optics.  The  simulations  randomly  generate  millions  of  x-ray  photons  from  a  particular  location  and  solid 
angle,  and  then  trace  them  through  the  optics  using  reflection  theory  and  classic  mechanics.  The  neces¬ 
sary  confidence  in  the  simulations  is  being  built  by  comparing  the  results  of  the  modeling  program  with 
experimental  data. 


For  single  fibers,  the  computational  speed  is  greatly  enhanced  by  a  reduction  to  two  dimensions  by  pro¬ 
jecting  the  trajectory  onto  the  local  fiber  cross-section.8  Reflections  are  computed  from  standard  tables.9 
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1.E  Purpose 

The  purpose  of  this  project  is  to  contribute  to  an  ongoing  effort  to 
use  the  new  innovative  technologies  in  X-ray  optics  to  develop  a 
digital  mammographic  imaging  system  with  high  resolution  and 
contrast  by  extending  the  capabilities  of  the  computer  simulation. 

Enhanced  computer  simulations  are  critical  because  manufactur¬ 
ing  large  polycapillary  optics  is  difficult,  expensive  and  time- 
consuming,  and  simulation  is  a  good  design  tool  to  predict  the 
performance  of  the  optics  and  make  it  clear  how  all  the  different 
parameters  work  together  and  impact  the  resultant  transmission, 
resolution  and  contrast. 

1. F  Technical  objectives 

The  specific  aims  of  this  traineeship  are: 

•  to  simulate  and  measure  trial  optic  setups,  to  compare  the 
simulated  results  with  experimental  data  and  to  determine  the 
best  design  for  building  the  best  direct  digital  mammographic 
system; 

•  to  further  the  author’s  knowledge  about  the  mammographic  systems  and  to  promote  a  future  career 
in  this  field. 

2.  Body:  Methods  and  Results 

2. A  Introduction  to  profile  and  surface  defect  corrections 

Real  glass  surfaces  are  not  perfect.  In  addition  to  absorption,  there  is  roughness,  which  has  short  wave¬ 
lengths,  of  the  order  of  microns  or  less,  waviness,  which  has  longer  wavelength,  and  bending,  which  has 
wavelength  longer  than  the  length  of  the  capillary.  Bending  is  described  by  the  radius  of  curvature  of 
the  fiber,  R  in  Table  1.  Waviness  results  in  the  channel  wall  being  tilted  with  respect  to  the  axis  of  the 
fiber,  and  therefore  changes  the  angle  of  incidence.  It  is  described  in  terms  of  a  distribution  of  tilt  an¬ 
gles.  Surface  roughness  produces  both  diffuse  scattering,  for  which  the  reflected  angle  can  be  either 
larger  or  smaller  than  the  corresponding  incident  angle,  and  transmission  into  the  glass.  The  roughness 
correction10  used  in  previous  simulations"  is  described  in  terms  of  a  roughness  height,  z  and  a  correla¬ 
tion  length,  s. 

2.B  Current  progress 

To  provide  a  framework  for  assessing  the  results  of  the  first  year  of  the  project,  the  original  statement  of 
work  for  the  full  three-year  project  is  reproduced  below. 

Proposed.  Statement  Of  Work 

Task  1  Developing  and  testing  the  simulation  software,  (months  1-18) 

Task  2:  Perform  experiment  and  geometric  simulations  on  pre-patient  collimating  optic 
at  a  wide  range  of  energies  (month  19-23) 

2. 1  Measure  transmission  as  a  function  of  x-ray  energy. 
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Figure  3.  The  generation  of  a  ran¬ 
dom  point  (of  a  photon)  in  a  circle. 
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2.2  Measure  transmission  uniformity. 

2.3  Measure  dependence  of  transmission  on  x-ray  entrance  angle. 

2.4  Perform  geometric  simulation  on  transmission  spectra. 

Task  3:  Perform  experiments  on  collimating  optic  at  mammographic  energies. 

(month  24-26) 

3.1  Measure  exit  angle  divergence. 

3.2  Measure  contrast  enhancement. 

Task  4:  Perform  experiments  and  geometric  simulations  on  post-patient  anti  scatter 
optic,  (month  27-30) 

4.1  Measure  central  part  and  whole  optics  as  a  function  of  photon  energy 

4.2  Measure  transmission  uniformity 

4.3  Measure  dependence  of  transmission  on  x-ray  entrance  angle 

4.4  Perform  geometric  simulation  on  transmission  spectra 

Task  5:  Perform  experiments  on  post-  patient  anti  scatter  optic  at  mammographic 
energies  (month  3 1  -33) 

5.1  Measure  transmission  energy  as  a  function  of  photon  energy 

5.2  Measure  transmission  uniformity 

5.3  Measure  transmission  entrance  angle  dependence 

5.4  Measure  scatter  fraction 

5.5  Measure  contrast  enhancement 

5.6  Measure  uniformity  of  magnification  factor 
Task  6.  Design  assessment  (month  34-36) 

In  the  first  year  of  this  project,  we  developed  and  tested  a  more*  physical  model  and  compared  the  simu¬ 
lation  results  for  a  variety  of  straight  capillaries  to  measured  data  in  the  energy  range  from  10  to  80  keV. 
Using  this  analysis,  the  number  of  required  parameters  in  this  energy  range,  including  20  keV  for  mam¬ 
mography,  is  reduced  from  four  to  two,  bending  curvature  and  waviness.  It  was  not  found  to  be  neces¬ 
sary  to  include  roughness  effects.  Although,  it  took  more  than  a  year  to  develop  this  unexpected  new 
model  and  perform  its  implementation  for  single  fibers  and  lenses,  it  was  worthwhile  because  we  now 
have  a  more  physical  simulation  model  which  increases  the  time  efficiency  for  future  modeling. 

Sections  2C-2E  are  the  comparison  of  the  new  model  with  the  previous  model  and  the  comparison  of 
their  simulation  results  with  the  experimental  data  for  borosilicate  glass  optics.  Section  2F-2G  presents 
new  data  and  simulations  applied  during  the  second  year  to  lead  glass  optics.  During  the  second  year  the 
principal  investigation  was  changed  from  Hui  Wang  to  Cari. 

2.  C  Improvement  in  generating  x-ray  photons 

As  shown  in  Figure  3,  to  produce  a  random  point  P  in  a  circle,  we  can  randomly  generate  either  (r,  0)  or 
(x,  y).  Previous  simulations' 1  used  uniformly  generated  0  in  the  range  (0,  2n)  and  r  in  the  range  (0,  R), 
where  R  is  the  radius  of  the  x-ray  source,  capillary  channel  or  fiber  input  area.  This  work  uniformly 
generated  x  and  y  coordinates,  with  the  limitation  x2  +  y2  <  R2,  which  is  closer  to  an  isotropic  distribu¬ 
tion. 
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Improvements  in  waviness  correction  model 


Waviness  of  a  capillary  surface  is  oscilla¬ 
tions  with  wavelengths  shorter  than  the 
capillary  length  and  longer  than  the 
wavelength  of  the  roughness.  The  wavi- 
ness  effect  is  modeled  as  a  random  tilt  of 
the  glass  wall,  which  results  in  a  random 
change  in  the  reflection  angle  on  each 
bounce.  If  the  random  tilt  angle  is  <j)  and 
the  incident  angle  is  0,  then  the  reflection 
angle  with  respect  to  the  nominal  surface 
will  be  9+2<|>,  as  shown  in  Figure  4.  Pre¬ 
vious  simulations  used  a  uniform  distri¬ 
bution  model;  this  work  uses  a  tilt- 
corrected  normal  distribution  to  generate  <j). 


Figure  4.  Scheme  of  an  x-ray  beam  interactions  with  a 
randomly  tilted  surface. 


2.D.1  Uniform  distribution  model 

Previous  simulations,  referred  here  as  Model 
Ml ,  assumed  a  uniform  distribution  of  2(j) 
from  — A0ijig\  to  A0niax  if  0  ^  AO ninx ■  or  from  -0 
to  A0max  if  0  <  A9max.  The  maximum  random 
tilt  angle  A9max  is  the  adjustable  parameter 
which  describes  the  amount  of  waviness  of  the 
capillary. 

2.D.2  Tilt-corrected  normal  distribution 
model 

2.D.2.1  Normal  distribution 


Figure  5.  Three  surfaces ,  OA/,  OA2  and  OA3,  with 
different  till  angles,  <f>\,  fc  and  fc,  respectively,  from 
the  nominal  surface  OO Here,  -6  < </>j<0, 


For  this  work,  named  as  Model  M2,  it  is  as-  ff  ™  ^ ™  ^ 

sumed  that  these  tilt  angles,  <j>,  are  normally  1  cos  ^ 1  2  cos  ^2  cos  ^  3  ~  ' 

distributed  in  the  range  (-n/2,  n/2).  with  the 

mean  value  equal  to  zero.  For  high  quality  optics,  the  standard  deviation  of  this  normal  distribution,  a,  is 
much  smaller  than  the  critical  angle,  0C.  The  probability  distribution  of  tilt  angles,  <J>,  is 


•  -s fin 


2.D.2.2  X-ray  impact  tilt  correction 

In  this  work,  Model  M2,  consideration  was  taken  of  the  fact  that  the  surface  tilt  angle  will  affect  the 
probability  of  x-ray  impact  on  that  surface.  Taken  to  extremes,  a  surface  region  perpendicular  to  the 
beam  is  much  more  likely  to  intercept  the  beam  than  a  surface  region  parallel  to  the  beam.  Figure  5  dis¬ 
plays  three  surfaces,  OAi,  OA2  and  OA3,  with  different  tilt  angles,  <|>i,  <|>2  and  fo,  respectively,  from  the 
nominal  surface  OO’.  The  projections  onto  the  nominal  surface  for  the  three  surfaces  are  equal, 
OA|*cos(|>i  =  OA2*cos<|)2  =  OA3*cos<(>3  =  OO’.  The  probability  of  incidence  of  a  parallel  x-ray  beam  with 
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incident  angle  0  (with  respect  to  the  nominal  surface  00’)  hitting  these  tilted  surfaces  is  given  by  their 
corresponding  perpendicular  length  AjBj  (for  j  =  1,  2  or  3) 

P.  x  A,Bf  -OAf  -  sin(6^  +  <pf )=-——■  sin(0  +  <fij ), 
cos^. 

We  call  this  the  tilt-corrected  probability  distribution.  The  complete  description  of  this  distribution  is 

'Fsin<g  +  *>,  -<?<«*<* 

//(*)  =  cos^>  2 

o,  --<  (j)  <  -e  ^  ' 

L  2 

? 

where  0  is  the  incident  angle,  <(>  is  the  tilt  angle  and  F  is  a  normalization  constant. 

2.D.2.3  Tilt-corrected  normal  distribution 

Combining  the  normal  distribution  G(<)>)  with  the  tilt  correction  H(4>)  gives  (j)  for  a  certain  incident  angle 
0  as  follows: 

r  -dP 

„  yj  sin (@+</>)  n  ,  n 

Ke2a - - — -Q  <</)<—  (r, 

J{(/>)=\  cos^  2,  (6) 

o, 

where  0  is  the  incident  angle,  <|)  is  the  tilt  angle  and  K  is  a  normalization  constant. 

Noting  that 

2^=sin*+corf-ta* 

cos^  (7) 

and  that  both  the  incident  angle,  0,  and  the  tilt  angle,  <|>,  are  very  much  less  than  40  mrad,  we  use  the  ap¬ 
proximation 

it 1 

ce1*1  *{e  +  <f>),  -e  <</><- 

P{<j>)  =  \  2,  (8) 

0,  -  —  <d><-d 

l  2 

where  the  normalization  constant ,  C,  is 


9a  +  2 a 2 


to  speed  the  calculation.  P(<j>)  is  called  the  tilt-corrected  normal  distribution.  In  model  M2,  a  is  the  vari¬ 
able  that  describes  the  amount  of  waviness.  Two  examples  of  J(<j>)  and  P(<|))  are  shown  in  Figure  6  and 
Figure  7. 
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Fiber  Description 

Model  Ml 

Model  M2 

Fiber 

# 

Type 

Outer 

Diame¬ 

ter 

mm 

Channel 

Size 

pm 

Open 

Area 

Length 

mm 

z 

nm 

s 

(am 

R 

m 

A0ninx 

mrad 

R 

m 

a 

mrad 

A 

1 

0.5 

12 

65% 

105 

0.7 

6 

105 

0.4 

60 

0.225 

C 

3 

0.75 

22 

50% 

136 

0.5 

6 

125 

0.35 

225 

0.2 

D 

4 

4 

12 

55% 

130 

0.8 

6 

110 

0.285 

139 

0.125 

E 

5 

0.3 

4-5 

55% 

105 

0.7 

6 

28 

0.2 

31 

0.09 

F 

4 

4 

12 

55% 

130 

0.8 

6 

90 

0.45 

90 

0.18 

Table  1  Parameters  for  best-fit  simulations.  R  is  the  bending  radius.  For  Ml,  AOmaxis  the  amount 
of  waviness,  z  is  the  roughness  height  and  s  is  the  roughness  correlation  length.  For  M2,  a  is  the 
standard  deviation  of  the  waviness. 


Figure  6.  The  probability  distribution,  J( </>),  and 
its  approximation,  P (</>),  versus  tilt  angle  (/>  when 
incident  angle  9=0. 0001  rad,  standard  deviation 
<j  =  0. 0004  rad.  <fi  is  in  rad. 


Figure  7.  J(<j>)  and  P(f)  versus  lilt  an¬ 
gle  <f>  when  incident  angle  6  =  0. 009  rad,  stan¬ 
dard  deviation  <7=  0.0004  rad.  <f>  is  in  rad. 


2.E  Results 

Transmission  of  a  number  of  different  fibers  measured  previously1 1  in  the  energy  region  1 0-80  keV  are 
shown  in  figures  8-18.  The  descriptions  of  these  polycapillary  fibers  and  their  simulation  parameters  are 
shown  in  Table  1. 
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Figure  8  Transmission  spectra  Figure  9  Simulations  of  Figure  10  Simulations  of  trans- 

for  fiber  C  simulated  with  transmission  spectra  for  fiber  mission  spectra  with  waviness 
different  bending  curvature  C  with  only  waviness  com-  correction  using  different  numbers 
alone  using  models  Ml  and  pared  with  the  experimental  of  photons,  compared  with  the  ex- 
M2,  compared  with  expert-  data.  This  figure  shows  the  perimenlal  data  of  fiber  C.  The 
mental  data.  This  figure  shows  effects  of  waviness.  It  also  figure  shows  that  the  simulation 
that  the  simulations  with  shows  simulations  using  curves  using  10000 photons  are 

bending  alone  can  not  fit  the  waviness  alone  do  not  fit  the  much  smoother  than  those  using 
data,  and  indicates  that  data.  The  simulations  do  not  1000 photons  because  of  the  re¬ 

bending  is  not  the  only  factor  include  the  roughness  or  *•  duction  in  statistical  fluctuation, 
to  cause  the  transmission  bending  correction, 

drop. 
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The  experimental  data  and  simulation  results  from  the 
two  simulations  for  a  type  3  polycapillary,  Fiber  C,  are 
compared  first  to  demonstrate  the  fitting  processes,  to 
interpret  the  waviness  and  bending  effects,  and  to  un¬ 
derstand  the  final  results. 

2.E.1  Bending  effects 

A  slight  bending  can  dramatically  reduce  the  transmis¬ 
sion  of  high-energy  photons  because  of  the  small  criti¬ 
cal  angle  at  high  energies.  Figure  8  shows  a  compari¬ 
son  between  experimental  data  and  simulations  with 
different  bend  radii.  The  difference  in  the  bending  ef¬ 
fect  between  Ml  and  M2  is  due  to  the  different  algo¬ 
rithms  used  in  generating  x-ray  photons  as  discussed  in 
section  2.C.  Ml  assumes  more  photons  hit  the  center  of 
the  fiber  as  well  as  the  center  of  each  channel.  M2 
generates  photons  uniformly  across  the  channel,  fiber 
and  source. 

For  both  models,  a  bend  of  R  =  80  m  reduces  the 
transmission  too  much  at  high  energies.  Since  any 
waviness  or  roughness  correction,  like  the  bending  cor¬ 
rection,  will  reduce  the  transmission,  R  must  be  larger 
than  80  meters.  We  set  Rmjn  =  80  m,  which  is  useful  in 
determining  the  best-fit  parameters  for  the  model  M2. 


Photon  energy(KeV) 

Figure  11.  Simulated  transmission  spectra 
using  model  M2  with  fixed  waviness  (0=0.2 
mrad)  and  different  bending,  compared  with 
the  experimental  data.  The  simulation  curves 
shown  are  just  a  few  representative  cases  in 
the  bisection  process.  For  M2,  the  best-fit 
parameters  are  <7=  0.2  mrad  and  R  =  225 
m;  for  Ml,  they  are  z  =  0.5  nm,  s  =  6  pm, 

A  0max  =  V.35  mrad  and  R  =  125  m. 


Source  position(mm) 


Figure  12.  Simulations  of  source  scan  curve  us¬ 
ing  model  M2  with  different  waviness  (o)  com¬ 
pared  with  the  experimental  data.  The  larger  the 
waviness,  the  narrower  the  simulated  scan  curve. 
The  simulation  curves  with  the  best  fitting  pa¬ 
rameters  of  Ml  and  M2  obtained  in  Figure  11 
are  also  shoM’n  in  this  figure  for  reference. 


Source  position(mm) 

Figure  13  Simulations  of  source  scan  curves 
using  model  M2  with  no  waviness  and  differ¬ 
ent  bending  compared  with  the  experimental 
data  at  10  keV.  The  simxdated  scan  curve 
does  not  become  narrower  even  when  the 
bending  radius  R  =  20  m.  The  word  "ideal  ” 
refers  to  a  simulation  with  R=ooand  no 
waviness. 
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Clearly,  bending  alone  cannot  explain  the  experimental 
data,  especially  at  low  energies. 

2.E.2  Waviness  effects 

To  explain  the  data  at  energies  from  20  to  60  keV,  an¬ 
other  correction  is  needed.  Roughness  does  not  have  a 
major  effect  in  this  energy  range1 '.  Waviness  is  ex¬ 
pected  to  be  the  major  factor  in  this  energy  range. 

Figure  9  shows  that  the  transmission  becomes  lower 
for  all  the  energies  and  the  high-energy  transmission 
drops  faster  with  increasing  a  gets  larger.  It  also  shows 
that  waviness  alone,  like  bending  alone,  can  not  ex¬ 
plain  the  experimental  data.  In  the  simulation,  a  =  0.3 
mrad  and  0.4  mrad  are  definitely  over-correcting  at  low 
energies,  so  we  can  set  0.3  mrad  as  the  maximum  of  a 
(crniax)  for  Fiber  C.  The  waviness  parameter  A0  for 
model  Ml  is  larger  than  the  parameter  a  for  model  2 
because  the  tilt  correction  in  model  M2  makes  it  more 
probable  that  the  x-ray  will  intersect  the  high  angle  sur¬ 
faces. 

The  fluctuation  of  the  simulation  curve  in  Figure  9 
when  a=0.25  mrad  is  mainly  due  to  the  relatively  small 
number  of  photons,  1000,  used  in  the  simulation.  Figure 
photons  changes  from  1000  to  10000. 


Source  position(mm) 

Figure  14  Simulations  of  source  scan  curves 
using  model  M2  with  fixed  waviness  (a  =  0.2 
mrad)  and  different  bending  compared  with 
the  experimental  data.  The  simulated  scan 
curve  does  not  apparently  change  even  when 
the  bending  radius  R  =  20  m.  With  and  without 
the  bending,  simulations  using  a  =0.2  mrad  fit 
the  experimental  data  quite  well. 


10  shows  the  difference  when  the  number  of 


Using  a  bisection  algorithm  with  Rmin  and  amax,  R  and  a  were  applied  together  and  compared  to  the  data 
to  get  the  best  fit  parameters,  as  shown  in  Figure  1 1 .  The  best  fit  parameters  for  fiber  C  using  model  M2 
are  R  =  225  m  and  ct  =  0.2  mrad. 

2.E.3  Verification  using  the  source  scan  curves 

Figure  12  shows  the  source  scan  curve,  that  is  transmission  as  a  function  of  lateral  source  displacement, 
and  simulation  results,  for  increasing  waviness  for  Fiber  C  at  10  keV.  The  simulated  source  scan  slowly 
becomes  narrower  when  a  gets  larger,  however  that  change  is  too  small  to  be  used  to  determine  a. 
Figure  13  displays  simulations  for  Fiber  C  at  10  keV  with  different  R  and  no  waviness.  Figure  14 
shows  the  effect  of  different  R  with  a  fixed  a.  From  these  two  figures,  it  is  apparent  that  changing  R 
from  20  meters  to  infinity  does  not  have  any  visible  effect  on  the  fitting  curve.  Therefore,  we  can  not  use 
the  simulations  of  the  source  scan  curves  to  determine  either  a  or  R,  but  they  can  be  used  to  verify  a. 

For  Fiber  C,  as  shown  for  10  keV  in  Figure  12,  at  higher  energies  in  Figure  15,  and  as  a  function  of  en¬ 
ergy  in  Figure  1 1,  the  simulation  curve  of  a=0.2  mrad  fits  the  experiment  data  very  well  with  one  fewer 
fitting  parameter  than  the  best  fit  of  Ml. 
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20KeV  40KeV:  59KeV:  80KeV: 

*  data  *  data  *  data  *  data 

Ml:  Best  fit  Ml:  Best  fit  Ml:  Best  fit  —■—Ml:  Best  fit 

•  •  O  -M2:ct=0.2  ~ -o ~  M2: ct=0.2  -- O--M2:o=0.2  - -O- -  M2:a=0.2 

— A--M2:a=0.2,R=225  A  ’  M2:<j=0.2,R=225  ^  M2:o=0.2,R=225  '  A  •  M2:rr=0.2,R=225 


Source  position(mm)  Source  position(mm)  Source  position(mm)  Source  position(mm) 

Figure  15  Simulated  source  scan  curves  using  Ml  and  M2  compared  with  experimental  data  at  four 
different  photon  energies.  They  all  fit  the  data  quite  well  For  M2,  this  figure  also  verifies  that  the 
bending  has  no  visible  effect  on  the  source  scan.  The  best-fit  parameters  for  Ml  are  R  =  125  m,  A6max 
=  0.35  mrad,  s  -  6  pm  and  z  =  0.5  nm. 

2.E.4  Best-fits  of  all  the  fibers  and  analysis  of  the  results 


Using  the  algorithm  described  above,  the 
best-fit  parameters  of  M2  for  Fiber  A,  D 
and  E  were  obtained.  The  best-fit  parame¬ 
ters  for  both  Ml  and  M2  are  shown  in  Ta¬ 
ble  1.  The  best- fit  curves  of  Ml  and  M2 
along  with  the  experimental  data  are  shown 
in  Figures  16-18.  The  fibers  are  described 
in  Table  1. 

In  Figure  1 7,  the  transmission  for  fiber  A 
shows  a  rapid  drop  for  energies  above  30 
keV.  Although  fiber  D  has  lower  fractional 
open  area  than  fiber  A,  its  transmission  ex¬ 
ceeds  that  of  fiber  A  at  energies  above  30 
keV.  This  is  because  fiber  A  is  thin  (0.5 
mm  in  outer  diameter)  and  flexible,  there¬ 
fore  difficult  to  keep  straight  in  the  meas¬ 
urement  apparatus.  Model  M2,  which  can 
vary  only  bending  and  waviness,  requires  a 
much  sharper  bend  for  fiber  A  than  for  fi- 


7Q  Fiber  C:  ^  data,-  ••-'Ml:  best  Fit; A  •  M2:cr=0.2,R=225 
Fiber  D:  X  data,' Ml:  best  fit, ’O'"  M2:<r=0. 125,  R=139 
Fiber  E— +—  data,-  -T~  M 1 :  Best  fitr<>—  M2:a=0.09,  R=3 1 


o-l - - - , - - - , - . - , - - - , - . 

0  20  40  60  80 

Photon  energy  (KeV) 

Figure  16.  Simulations  of  transmission  spectra  of fiber 
C,  fiber  D  and  fiber  E  with  their  best-fit  parameters 
compared  with  the  experimental  data.  All  parameters 
related  to  this  figure  are  listed  in  Table  1. 
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ber  D.  This  bending  is  more  significant  at  high  energies,  where  the  critical  angles  are  smaller. 

Transmission  curves  for  fiber  C  and  D  are  similar  in  shape  and  are  nearly  flat  up  to  60  keV,  as  shown  in 
figure  15.  However,  the  smaller  channel  size  for  fiber  D  results  in  the  simulation  being  less  sensitive  to 
bending.  Therefore  a  larger  bending  curvature  (smaller  radius)  is  required  for  fiber  D.  The  waviness 
correction  for  fiber  D  must  therefore  be  smaller  to  keep  the  transmission  about  the  same  as  for  fiber  C. 

Fiber  E  is  the  thinnest  fiber.  Table  1  shows  that  it  has  the  largest  bending  correction  because  of  its  flexi¬ 
bility.  Its  transmission  curve  is  flat  up  to  40  keV.  This  is  due  to  its  very  small  channel  size,  as  discussed 
in  the  comparison  of  fiber  C  and  fiber  D.  However,  if  the  channel  size  is  too  small,  as  in  fiber  E,  it  also 
results  in  more  reflections  being  needed  for  a  photon  to  traverse  the  fiber  and  may  have  introduced  other 
defects  such  as  blocked  channels.  This  is  why  the  transmission  is  only  40%  for  the  energies  below  40 
keV  although  the  open  area  is  around  55%. 


The  high  transmission  and  the  simulation  results 
good.  The  bending  radius  is  above  130  meters 
for  type  3  and  type  4  capillaries.  It  is  hopeful 
that  we  can  further  improve  the  high-energy 
transmission  performance  of  polycapillary  fi¬ 
bers  by  decreasing  to  an  optimum  channel  size, 
and  making  them  more  rigid. 

Because  Model  M2  is  physically  plausible, 
gives  reasonable  results,  and  matches  the  ex¬ 
perimental  data  very  well,  as  shown  in  figures 
16  and  17,  we  can  use  it  to  simulate  other  poly¬ 
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Figure  18.  Simulations  of  transmission  spectra 
of  fiber  D  (unexposed)  and fiber  F  (exposed) 
with  their  best-fit  parameters  compared  with 
the  experimental  data.  Layer  =  35  means  that 
there  is  a  35  micron  glass  layer  added.  All  pa¬ 
rameters  related  to  this  figure  are  listed  in  Ta¬ 
ble  1  except  the  glass  layer. 


show  that  the  quality  of  the  capillary  fibers  is  quite 
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Figure  1 7.  Simulations  of  transmission  spectra  of 
fiber  A  and  fiber  D  with  their  best-fit  parameters 
compared  with  the  experimental  data.  All  parame¬ 
ters  related  to  this  figure  are  listed  in  Table  1. 


capillary  optics.  Compared  to  Ml,  which  uses  four 
parameters,  M2  uses  only  two  parameters,  which 
makes  simulating  the  experimental  results  more 
straightforward,  easier  and  timesaving.  In  addition, 
roughness  does  not  have  much  effect  in  the  energy 
range  10-80  keV,  so  M2  gives  us  more  confidence 
with  fewer  parameters. 

2.E.5  An  analysis  of  a  radiation-damaged  fiber 

The  measured  transmission  spectra  and  source  scans 
provide  sensitive  tools  for  analyzing  capillary  quality 
This  method  was  used  to  analyze  the  effect  of  radia¬ 
tion  damage  on  fiber  F  that  has  the  same  type  and 
same  length  as  fiber  D.  Radiation  damage  has  been 
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observed  after  a  very  large  dose  of  x-ray  exposure.12  It  was  found  that  prolonged  exposure  to  synchro¬ 
tron  white  beam  radiation  could  cause  a  measurable  bend  for  a  thin  fiber.  This  bending  can  be  reduced 
by  annealing  and  by  holding  the  fiber  rigid  during  exposure.  Fiber  F  had  1 .8  MJ/cm2  exposure  before  the 
current  measurements.  There  is  no  visually  observable  bending  or  other  visible  damage. 

The  transmission  spectrum  and  its  simulation  for  fiber  F  are  compared  with  an  unexposed  capillary  of 
the  same  type,  fiber  D,  in  Figure  16.  The  measured  data  shows  a  significant  transmission  drop  above  30 
keV  and  below  10  keV.  The  simulation  indicates  that  the  exposed  capillary  has  more  bending  and 
waviness.  The  bending  radius  was  changed  from  139  m  to  90  m  and  the  waviness  ct  changed  from  0.125 
mrad  to  0. 1 8  mrad.  Although  these  profile  changes  have  a  significant  effect  the  transmission,  they  are 
still  very  small  and  not  visible.  The  transmission  drop  below  20  keV  can  be  explained  by  adding  35  fj.m 
x-ray-absorbing  glass  layer  as  shown  in  Figure  18. 13 

2.F  Measurement  and  analysis  of  leaded  glass  fibers 

As  discussed  in  the  previous  section,  the  geometric  simulation  model  M2  has  been  used  to  analyze  the 
experimental  data  for  number  of  borosilicate  single  capillaries.  The  analysis  shows  that  the  simulation 
results  fit  fairly  will  with  the  experimental  data. 

In  this  section,  we  will  analyze  the  experimental  results  of  a  number  of  leaded  glass 
single  capillaries  using  model  M2.  Leaded  glass  polycapillary  optics  could  be  a  candidate  for  scatter  re¬ 
jection  because  of  lead’s  superior  absorption.  Lead  glass  optics  could  be  shorter  and  therefore  easier  to 
manufacture.  Lead  glass  optics  may  be  the  fast  way  to  be  bring  capillary  optics  to  clinical  use. 
Transmission  of  a  leaded  glass  polycapillary  in  the  9-80  keV  energy  range  is  shown  in  figure  19.  The 
description  of  these  capillaries  and  the  simulation  parameters  listed  in  table  II.  The  characteristic  of  the 
transmission  spectra  leaded  glass  fiber  are  similar  to  the  radiation-damaged  borosilicate  fiber,  i.e,  the 
transmission  drops  off  at  both  lower  energies  and  higher  energies,  as  shown  in  figure  19.  In  addition  to 
the  waviness  and  bending  simulation  parameters,  the  simulation  has  to  include  a  leaded  glass  layer  of 
50  pm  thick  or  less,  to  fit  the  experimental  data. 


Table  II.  Fiber  description 


Model 

M2 

Model 

M2 

Fi¬ 

ber 

# 

Outer 

di¬ 

ameter 

Chan 

nel 

size 

m 

Open 

area 

% 

Length 

mm 

s 

jam 

R 

m 

W 

mrad 

Z 

nm 

Layer 

pm 

R 

m 

W 

mrad 

Layer 

pm 

A 

0.54 

12 

60 

60 

6 

33 

0.5 

2 

50 

42.5 

0.3 

54 

B1 

0.51 

12 

60 

60 

6 

52.5 

0.3 

2 

30 

59 

0.20 

5 

26 

B2 

0.51 

12 

60 

30 

6 

30 

0.3 

2 

35 

30 

0.20 

5 

37.5 

C 

0.54 

11 

50 

60 

6 

28 

0.5 

2.25 

30 

27.5 

0.26 

37 
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Figure  19.  Simulated  transmission  spectra  us¬ 
ing  model  Ml  and  M2  for  leaded  glass  single 
fiber  type  B1  of  length  60  mm  compared  to  the 
experimental  data. 


a  type  A  exp 
•  typeB2:  exp 


Figure  20.  Simulated  transmission  spectra  using  model  M2 
for  leaded  glass  single  fiber  A,  fiber  B  and  fiber  C.  All  pa¬ 
rameters  compared  with  experiment  data  and  related  to  this 
figure  are  listed  in  table  II. 


Figure  19  shows  that  both  simulations  model  Mi  and  model  M2j  fit  fairly  well  with  the  transmission 


■  type  B1:  exp9keV 

model  M2 


Source  Position(mm) 


Figure  21.  Simulations  of  source  scan  curve  using 
model  Ml  and  M2  compared  with  experiment  data 
for  leaded  glass  for  single  fiber. 


■  exp30kev 


SoLnoe  Rsition(rrm) 

Figure  22.  Simulation  of  source  scan  curves  using 
model  Ml  and  M2  compared  with  experimental  data 
for  leaded  glass  for  type  Bl. 
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■  exp  30  kev  spectrum  for  the  type  B 1  fiber 
•  exp  40  kev  once  a  layer  of  leaded  glass  is 

1-0  ‘  ^  -  A  exP  69  kev  added..  Since  model  M2  does  not 

g  2 ;ijl  T  exP  80  kev  include  the  surface  roughness 

’jg  0  8 "  *j  jj  \\  \  ■  model  M2  and  leaded  glass  fibers  have  very 

|  /  h  f  *!  model  M2  high  roughness,  it  is  expected 

ro  0  6  “  %/  jj  |\  f  model  M2  that  model  M2  does  not  fit  the 

£  /  /  / /  \\  \  \  . ~  model  M2  source  scan  data  well  at  energy 

|  04 '  1  /  J1  V\  \  «i  20  keV  as  well  as  model  Ml,  as 

|  /  •  //  \\  \  \  shown  in  figure  21.  Neither 

o  0  2 "  /m  /  )j  ▼\\  \  model  fits  the  data  well  at  9keV 

z  Vk  However,  model  M2  fits  the 

o.o  -  *»****#*$**+'  ▼♦♦♦•♦Im******  source  scan  data  fairly  well  at 

- r - - - 1 - ’ - 1 - 1  higher  energies,  as  shown  in  fi- 

figure  gure  22.  0  2  gure22. 

Source  Position(mm)  .. 

rinally,  the  measured  trans- 

Figure  23.  Simulation  of  source  scan  curves  using  model  M2  com -  mission  spectra  for  a  number  of 

pat  ed  with  experimental  data  at  four  energies  for  type  Bl .  leaded  glass  fibers  II  are  shown 

in  figure  20. 

The  models  fit  the  transmission  spectra  fairly  well,  with  the  parameters  listed  in  table  II. 


▼  A' 


figure  gure  22.  2  0  2 

Source  Position(mm) 

Figure  23.  Simulation  of  source  scan  curves  using  model  M2  com¬ 
pared  with  experimental  data  at  four  energies  for  type  Bl. 


2G.  Measurement  of  the  transmission,  output  divergence,  and  uniformity  of  collimating  lens. 

This  collimating  lens  is  3  cm  x  3  cm  on  the  parallel  end,  and  2.9  cm  x  2.9  cm  on  the  focusing  end.  The 
focal  length  of  the  lens  was  determined  by  measuring  the  transmission  at  various  distances  from  the 
source.  The  highest  transmission  determines  the  focal  distance,  which  was  found  to  be  193  cm  for  this 
lens,  as  shown  in  figure  27. 


E  35  keV 

•  E40keV 

*  E  45  keV 

▼  E  50  keV 
o  E55keV 
+  E60  keV 


190  191  192  193  194  195  196  197  198  199 
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Figure  27.  The  transmission  as  function  of 
source-lens  distance  at  various  energies  for 
new  collimating  lens  with  focal  distance  of 
193  cm. 


■-  S-L  1930mm 
•  S-L  1980mm 
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Figure  28.  The  measured  transmission  v-.v 
energy  at  193  cm  and  198  cm  source-lens 
distances. 


Study  of  Capillary  X-ray  Optics  for  Digital  Mammography  Wang,  Hui  and  Cari  20 

The  transmission  of  this  lens  was  measured  by  placing  the  source  at  focal  distance.  Figure  28  shows  that 
this  lens  has  transmission  close  to  45  %  at  20  keV. 

To  measure  the  output  divergence  of  this  lens,  a  high  quality  single  crystal  (100)  silicon  wafer  is 
placed  between  the  exit  end  and  the  detector.  The  output  divergence  of  this  lens  is  less  than  6  mrad  as 
shown  in  figure  29. 


1390  13.95  14.00  14.06  14.10  14.15  14  20  14.25  14.30  -30  -15  -10  -5  0  5  10  15 

V\fefer  angle,  degrees  r*9^  (mm)  Energy  45keV 


Figure  29.  Output  divergence  of  colli-  Figure  30.  Output  intensity  distributions  for 

mating  lens,  measured  by  rocking  a  sili-  new  collimating  lens  at  198  cm  source  -lens 

con  crystal  in  the  beam  distance. 

>  ’ 

The  uniformity  of  the  lens  was  measured  by  placing  an  aperture  in  front  of  the  lens  and  taking  aperture 
scans  with  the  source,  lens,  and  the  detector  are  fixed.  Surface  area  plots  of  the  intensity  as  a  function  of 
aperture  position  can  be  seen  in  figure  30. 

2.H  Measurement  of  a  prototype  post-patient  monolithic  optic 

Polycapillary  optics  can  be  manufactured  either  as  monolithic  optics,  shaped  from  a  single  large 
fiber,  or  as  multi  fiber  optics,  with  multiple  fibers  strung  through  grids.  Although  a  lot  of  progress  has 
been  made  on  manufacture  of  monolithic  optics,  there  are  still  technological  difficulties  in  producing 
good  quality  large  scale  monolithic  so  far.  However,  magnification,  with  resulting  enhancement  requires 
monolithic  optics. 

In  this  section  transmission  measurements  for  lens  J  as  prototype  monolithic  optic,  are  presented. 
The  basic  geometric  of  this  post  patient  monolithic  optic  were  measured  as  following: 

Length  of  the  optic  =129  mm 
Output  diameter=8.25  mm 
Input  diameter=6. 1  mm 
Channel  size=25  pm 
Focal  distance=500  mm 

The  x-ray  generator  used  in  the  experiment  is  a  low  current  microfocus  MS50  with  a  50  pm  spot  size, 
tungsten  target,  and  a  maximum  operating  voltage  of  100  kVp.  Because  of  its  high  background  scattered 
x-rays,  this  source  is  enclosed  in  6mm  thick  lead  box  and  therefore  the  source  can  not  be  moved  in  z  di¬ 
rection.  The  optic  to  the  source  distance  was  then  changed  by  moving  the  optic  in  the  z  direction. 
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Variation  of  acceptance  angle  of  the  optics  as  a  function  of  source  optic  distance  was  used  primarily  to 
determine  both  focal  distance  of  prototype  monolithic  optic  and  the  input  depth  field.  The  acceptance 
angle  is  the  full  width  of  the  source  scan  divided  by  the  source  to  optic  distance.  The  acceptance  angle  is 
a  minimum  at  the  focal  distance,  where  the  transmission  is  maximized.  Focal  distance  of  this  optic  was 
determined  by  moving  the  source  in  x  and  y  direction  to  maximize  the  transmission  at  certain  position  of 
the  optic. 

Figure  24  shows  that  the  focal  length  and  the  input  depth  filed  of  this  optic  are  500  mm  and  20  mm,  re¬ 
spectively.  This  measured  focal  length  is  in  good  agreement  with  focal  length  determined  using  the 
transmission  as  function  of  source  -optic  distance  curve  as  shown  in  figure  24.  Finally  the  transmission 
efficiency  of  this  fiber  was  measured  at  the  focal  point.  This  optic  is  fairly  good  since  its  transmission 
efficiency  is  20.8  %  at  20  keV,  18  %  at  23  keV  and  16.33  %  at  25  keV.  At  mammographic  energy  (20 
keV)  the  high  angle  transmission  of  this  optic  close  to  zero  as  shown  in  figure  26. 


Figure  24.  Measured  acceptance  angle  as  function  Figure  25.  Measured  transmission  vs  source  -optic 
of  source-optic  distance  for  optic  J  at  20  ke  V  distance  of  optic  J  at  20  keV. 
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Figure  26.  Measured  of  scan  source  for  optic  J  at 
focal  distance. 


3.  Conclusion 

A  new,  more  physical  model  for  describing  pro¬ 
file  defects  in  polycapillary  x-ray  optics  at  high 
energies  for  mammography  has  been  developed. 
The  simulation  model  M2,  a  geometric  simula¬ 
tion  with  waviness  and  bending  correction,  is 
physically  plausible  and  works  quite  well  in  ex¬ 
plaining  the  transmission  spectra.  It  uses  only 
two  parameters,  which  makes  simulating  the  ex¬ 
perimental  results  and  designing  lenses  straight¬ 
forward,  easy  and  timesaving.  According  to  the 
simulations,  waviness  and  bending  in  the  capil¬ 
lary  channel  profiles  can  be  particularly  harmful 
for  high-energy  photons  because  of  the  smaller 


critical  angle  at  high  energies.  Roughness  is  not  important  in  the  simulation.  It  is  expected  that  increas¬ 
ing  the  rigidity  and  straightness  of  the  fiber  and  reducing  the  channel  size  to  an  optimum  value  can  fur¬ 
ther  increase  the  transmission.  The  spectra  of  transmission  versus  energy  have  been  found  to  be  a  very 
sensitive  tool  in  capillary  quality  analysis,  and  transmission  versus  source  scan  curve  measurements  can 
also  provide  significant  information.  Lead  glass  shows  promise  in  allowing  for  shorter,  higher  transmis¬ 
sion  optics  which  will  easier  to  manufacture.  The  new  collimating  lens  with  long  focal  distance  of 
around  2  m  has  good  transmission  and  uniformity. 
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